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ON CALCULATIONS OF ZEROS OF VARIOUS L-FUNCTIONS 



By Hiroyuki Yoshida 



Introduction. As we have shown several years ago [Y2], zeros of L(s, A) and 
L^"^' (s, A) can be calculated quite efficiently by a certain experimental method. Here 
A denotes the cusp form of weight 12 with respect to SL{2, Z) and L(s, A) (resp. 
L*-^-'(s, A)) denotes the standard (resp. symmetric square) L-function attached to 
A. The purpose of this paper is to show that this method can be applied to a wide 
class of L-functions so that we can obtain precise numerical values of their zeros. ^ 

We organize this paper as follows. In §1, we shall describe basic features of 

our method of calculation, which is repeated applications of partial summation. 

In §2, we shall study the r-th symmetric power L-function L'-'^)(s, A) attached to 

^ ' A. Since the cases r = 1, 2 are discussed in [Y2], we shall exclusively treat the 

cn . cases r = 3, 4. In §3, we shall study the L-functions attached to modular forms 

of half integral weight. These L-functions do not have Euler products. Naturally 

the Riemann hypothesis fails for them; we shall find many zeros off the critical 

line, though major part of zeros lie on the critical line. We shall also calculate 

Q\ • the location of these zeros off the critical line. Though there is some hope to find 

^ . relations among zeros of L-functions of two modular forms which are in the Shimura 

^ '. correspondence, no explicit results came out so far. 

In §4, we shall study L-functions attached to Hecke characters of non-Ao type 
of real quadratic fields. D.A. Hejhal showed great interest to make experiments 
in this case, since coefficients are non-computable combinatorially; hence there 
is a slight possibility that the Riemann hypothesis may break down for these L- 
5^ ! functions. We have made experiments on 44 cases summarized in Table 4.3; so far 

no counterexamples are found. 

In §5, we shall study the Artin L-function attached to a 4-dimensional non- 
monomial representation of Gal (Q/Q). In §6, we shall discuss the controle of error 
estimates in our calculation. In §7, we shall consider the explicit formula for the 
L-function attached to a modular form of weight 8 with respect to ro(2). We shall 
compare both sides of the explicit formula numerically. In §8, we shall present 
sample programs to compute values of L-functions, which may be convenient for 
the reader. In §9, we shall formulate a conjecture which emerged during the process 
of our experiments. 



* During the final stage of writing this paper, the author was at MSRI supported in part by 
NSF grant PMS9022140. 

^ After the publication of [Y2], H. Ishii [Is] published a table of zeros of standard L-functions 
attached to modular forms for 15 cases. It also comes to the author's notice that a program of 
the calculation of zeros of L{s, A) is included in "Mathematica" package, following the method of 
[Y2]. 
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Most of sections have attached tables to show results explicitly. Concerning 
actual computations, we have used "UBASIC " created by Y. Kida. (It was not 
available when we wrote [Y2].) The calculation was done by personal computers 
which are not necessarily so fast. However our experiments extended over long 
time (about three years) and UBASIC is quite fast (compared with some other 
softwares) for numerical calculations, the author thinks that our tables are fairly 
extensive. 

A motivation in these calculations has been to find non-trivial "functorial" prope- 
rties which may exist among zeros of L-functions, as was hinted in [Yl]. Though 
our experiments are not successful in this regard, conjectures stemmed from them 
are formulated in §9. 

We can pursue these calculations still further. The topics which may be included 
in this paper are: 

1) The Hasse-Weil zeta functions of algebraic curves, for example y'^ = x^ — x + 1. 

na — [na] — 1/2 
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The Dirichlet series 2, ~ studied by Hecke [H], where ct is a real 



n=l 

irrational number. 

3) Applications to Riemann-Siegel type formulas. 

4) Calculations of critical values of L-functions. 

Our results on these topics are still fragmentary, so the full discussion should be 
postponed to future occasions. 

Notation. For a complex number z, we denote by ^{z) (resp. ^'(-z)) the real 
(resp. imaginary) part of z. The letter q stands for exp(27rV— 12;) when it is clear 
from the context. For modular forms, we follow the notation in Shimura [Shi]. 

§1. An overview on our method of calculations 
Let 



(1-1) Hs) = J2 



00 
n=l 



be a Dirichlet series which is absolutely convergent when K(s) > a for some a > 1. 
In this paper, we shall consider only such L{s) which can be analytically continued 
to the whole complex plane as an entire function and satisfies a functional equation 
of the form 



(1.2) R{k-s)=KR{s) 

Here k is a constant of absolute value 1, /c > 0, 



R{s)=N'Y[r{b,s + c,)L{s) 



i=l 



with N > 0, bi > 0, Ci e C. We note that (1.2) is equivalent to 

(1.3) R{k-s)=KR{s) 
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if a^ G R, Ci G R for all n and i. Put k ^ = k^ with some Ki G C. By (1.2), we 
have 

(1.4) KiR{s) G R if ^{s) = k/2. 

Take any 6 > 0. For T > 0, let N{d, T) denote the number of zeros of R{s) 
counted with multiplicity in the domain 

-6 < 3f?(s) <k + S, < %s) < T. 

Let D be the rectangle whose vertices are —5, k + d, k + d + iT, —d + iT and let C 
denote the contour dD taken in positive direction. By the argument principle, we 
have 

(1.5) N(5,T) = — [ ^Mds, 
^ ^ ^''2'KiJc R{s) 

assuming that neither zeros nor poles of R{s) lie on C. Let Ci denote the portion 
of C from k/2 to k/2 + iT. By the functional equation (1.2), we obtain 

m 

(1.6) N{6, T) = TT-^Aarg R{s) = 7r-^A(arg A^" JJ r{biS + c^))+ 7r-^A(arg L(s)), 

i=l 

where Aarg denotes the variation of the argument on Ci, i.e., from s = k/2 to 
k/2 + iT along k/2 to k + 5, k + 5 to k + 6 + iT, k + 5 + iT to k/2 + iT. Set^ 

m 

7?(T) = Aarg (iV^ J] r(6,s + q)). 

Assume 

3f?(Q) > -b^k/2 for 1 < i < m. 

Then since bi > 0, A^^ Ili^i r(6iS + q) has neither zeros nor poles in the domain 
K(s) > k/2. Hence '&{T) is equal to the variation of the argument of A^'^ YllLi ^{biS+ 
Ci)) on the line segment [k/2, k/2 + iT]. We note that ^^(T) can be computed in 
high precision very easily using Stirling's formula (cf. [WW], p. 252) combined 
with the relation r(s + 1) = sr(s). We obtain 

(1.7) N{S, T) = 7r-M(T) + 7r-iA(arg L{s)). 

Now let us consider the case when R{s) has zeros in (—5, k + 6). Let r denote 
the number of zeros of R{s) i.e., of L{s), in this interval counted with multiplicity. 
Then (1.7) holds with the modification 

(1.8) N{S, T)-^ = Ti-^^T) + 7r-iA(arg L{s)). 

Here A(arg L{s)) is counted by dividing Ci into a finite number of paths removing 
real zeros of L(s) and summing the variations of the argument of L(s) on each of 



■^When it is clear from the context, we shall use i?(T) for the "phase factor" of this type in the 
following sections without further explanation. 
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them. The vaUdity of (1.8) can be seen by modifying C by small semi-circles which 
detour the real zeros of L{s). 

Throughout the paper, to compute L{s), we shall employ our method given in 
[Y2], which is repeated applications of Abel's partial summation. Set 



^n '•^n: "U-n ^ 



and define Sn , Un recursively by 



N 

Put Sp^ = Yl Snuy. Then we have 

n=l 

(1.10) s<;'=sr"-4'"r+v. 

As we have seen in [Y2], in several cases, S]^ seems to approximate L{s) amazingly 
well when we choose A^ and / sufficienly large. In the succeeding sections, we shall 
present various types of L-functions which can be treated in more or less similar 
fashion. The efficacy of our method seems to depend strongly on the arithmetical 
nature of the coefficients a„ of a Dirichlet series L{s). 

We shall conclude this section by technical remarks concerning actual compu- 
tations of Sj^ . As the first step, we should construct a table of an- For Dirichlet 
series considered in this paper, this step can be achieved rather easily. Since we 

can compute SpJ from S}^ by (1.10), the computation of S}^ = ^ ann~^ is the 

n=l 

substantial and the most time consuming part of our calculation. However usu- 
ally Sn becomes very large and Un very small when / increases. Therefore it is 
indispensable to perform the actual computation in high precision. For Wn , the 
following formula (1.11) should preferably be used than to compute it directly from 
the definition. 



'v-E(E (-!)"■ (1) 

k=l m=l ^ ^ 



(1.11) k=l m=l ^"V A^>/ > 1. 

,s(s+^2iLii±i:ill/v-'t 

If we replace Y^=i by Yk=ii ^^^ error is less than '^^ ^ (-^IIxm^ \j^)^^^\^~^\ if 
3fJ(s) > -L-1. 



§2. L-functions attached to symmetric tensor representations of GL{2) 

Let f{z) = Xl^i '^nC^'^'"^ G Sk{SL{2,Z)) be a normalized common eigenfunc- 
tion of Hecke operators. The L-function L(s, /) = Xl^i Cn^~'^ attached to / 
converges absolutely when ^{s) > -^^ and has the Euler product 



L{sJ) = l[{l-c,p-^+p''-'- 



■2s\-l 



Put 

l-CpX+ p^-^X^ = (1 - apX){l - (3pX) 

with ctp, Pp G C, where X is an indeterminate. For a positive integer r, we define 
an Euler product 

(2.1) lW(., /) = n[(i - «pp-^)(i - cx;-'PpP-n ■ ■ ■ (1 - /3>-^)]-' 



which converges absolutely when ^{s) > ^ ^ +1- K is conjectured that L('^)(s, /) 
can be analytically continued to the whole complex plane as an entire function and 
satisfies a functional equation. The conjectural functional equation of L^'^'{sjf) 
takes the following form (cf. Serre [Se]). If r is odd, put r = 2m — 1, 

m—l 

(2.2) i?W(., /) = (27r)— ^ n r(« - ^(^ - ^))L^'^Hs, /), 

i=0 



(2.3) Cr = (v/Zi)-+(fc-i)-'. 
If r is even, put r = 2m, 

m—\ 

R^^\s, f) = 7r-^/2(27r)— ^( J] r(. - z(fc - 1))) 

(2.4) *=o 

(2.5) e,. = (^/Zi)"^+(fc-i)"^("^+i)+5^ 

where 6 = (resp. 1) if to is even (resp. odd). Then the functional equation 

(2.6) i?W(s, /) = e,i?M(r(fc -l) + l-sj) 

is predicted. A quick way to see (2.6) is as follows. Let Mf be the motive of rank 
2 over Q attached to /. We see that the Hodge realization of Mf corresponds to 
the two dimensional representation 

p - Ind(V'fc; Wc -^ M^R,c) 
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of Wji,c- Here Wc = C^, Wrc is the Weil group of C over R and ipk is the 
quasi-character Vfc(a^) = x~^''~'^^ of Wc- Let a^ : GL{2) — > GL{r + 1) be the 
symmetric tensor representation of degree r and put p^ = Cr o p. Then we find 



(2.7) p, ^ ©^-ilnd(a; ^ a;- 


-{r-^){k-l)--^(k-l). ^^ ^ W^c)^ T 


p,-©— ^Ind(x- 

©{(x — > X' 


_^^-(r-^)(fc-l)--^(fc-l). ^^ _^ ^^ ^ 

-m(fc-i)(gg^^^m(fc-i)^Q^|^ r = 2m, 



2m -1, 



where t denotes the transfer map from VFr^c to R^. The gamma factor and the 
constant e^. of the functional equation can be calculated as the usual gamma factor 
and the constant attached to the representation pr of VFr,c; hence we obtain (2.2) 
~ (2.6). 

We refer the reader to Shahidi [Shal], [Sha2] for what are known on these sym- 
metric power L-functions, in more general cases. 

Let A(2) = Qnr=i(l - (f'f^ e Si2{SL{2,Z)), q = e^'^v^^ The calculation 
of zeros of L*-^-'(s, A) for r = 1, 2 is given in [Y2]. We consider the case r > 3. 
To compute L^'^^s^ A), we modify our summation method slightly in the following 
way. Fix r, choose v = Vr > and set 



lW(s, A) = Y^^nfi-' = 5^(a,n-^; 



n-("-"l 



n=l n=l 

Put 

and define s^ , Un recursively by (1.9). We set 5"^ = ^„=i -Sn "n • It turns out 
that a suitable choice of v depending on r yields good results. We can interpret 
this as the neutralization of the effect of extremely large value of sh, and extremely 
small value of Wn . 

As the first example, let r = 3. We take v = 8. For s = 17 + it, t = 20, the 
values of 

Rj = K(exp(z^?(t))>5}^^), Ij = $5(exp(z^?(t))5{^^) 

are given in Table 2.1. In Table 2.2, we give the values of t„, the n-th zero of 
L(3)(s, A), s = 17 + zt on the critical line for < t < 40. 

Next we apply our summation method to L^'^\s,A) taking v = 12. For s = 
^ + it, t = 10, the values of 

Rj = ^{exp{i^{t))Sii^), Ij = '^{exp{i^{t))S'ji^) 

are given in Table 2.3. In Table 2.4, we give the values of Un the n-th zero of 
L^'^\s, A), s = ^ -|- i-u on the critical line for < w < 20. 

We can see, by the same technique as will be given in §3 and §4, that the Riemann 
hypothesis holds for L^^\s,A) (resp. L^^'(s,A)) in the range < Q'(s) < 40 (resp. 
< '^{s) < 20 ) and that the zeros 17 + itn (resp. ^ + iun) are simple. 

§3. Modular forms of half integral weight 
6 



Put 

e{z) = J2 e(^'^) = 1 + 2 ^ Q^', viz) = e{z/24) J] (1 " eM). 



oo 
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n£Z n=l n=l 

By Shimura, [Sh2], (4.1), we have 

dim58(ro(2)) = 1, dim59/2(ro(4)) = 1 

and {7]{z)7]{2z))^ (resp. e{z)-^7]{2zy^) spans S8{To{2)) (resp. S9/2{To{4:))). Put 

CXD OO 

/(z) = ivizU2z)f = J2 ^nq\ 9{z) = 9{z)-'v{2zr = J^ ^^q^ 

n=l n=l 

CXD CXD 

Hs, f) = Y^ ar^n~% L{s, g) = ^Cr, 



n 



n=l n=l 



R{s, /) = r/^{2n)-^T{s)L{s, /), i?(s, g) = r{2n)-^T{s)L{s, g). 

Then / and g are in the Shimura correspondence; L(s, /) and L(s, g) can be ana- 
lyticaUy continued to the whole complex plane as entire functions and satisfy the 
functional equations 

(3.1) R{s, f) = R{8 ~sj), R{s, g) = R{9/2 -s,g). 

This example is described in detail in [Sh2], §4. For t > 0, let ^^/(t) (resp. ^?g(t)) 
denote the variation of the argument of 2*/^(27r)~*r(s) (resp. 2*(27r)~*r(s)) from 
4 to 4 + zt (resp. 9/4 to 9/4 + it). 

For L(s, f), s = 4 + it, t = 100, the values of 

Rj = 3?(exp(z^?/(t))5{^)), Ij = '^{eMi^f{t))S^^^ 



N 



are given in Table 3.1. 
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For L(s, g)j s = J + it, t = 100, the values of 



Rj = 3f?(exp(z^?,(t))5}^)), /, = Q{exp{i^g{t))S'}^^) 

are given in Table 3.2. 

By our method, we can compute zeros of L(s, /) and of L(s, g) on the critical line 
with sufficient accuracy observing sign changes of e*''-'"*^*-'L(4+zt, /) and 6^"^^^^' L{j + 
it,g). In Table 3.3, we list the n-th zero tn of L{s,f), s = 4 + zt in the range 
< t < 100. In Table 3.4, we list the n-th zero Un of L{s,g), s = | + i-u in the 
range < w < 100. 

Now let us examine the Riemann hypothesis for L(s, /). We see f{iy) > for 
y > by the product expansion of the ry-function. By the integral representation 

/>CXD 

(27r)-T(.)L(.,/)= / f{zy)y^-Uy, 

Jo 

7 



we see that L(s, /) > for s > 0. For T > 0, let N(T) denote the number of zeros 
of L(s, /) counted with multiphcity in the domain 

|3f?(s)-4| < 1/2, 0<S(s)<T. 

By (1.7) taking 5 = 1/2, we have 

(3.2) N{T) = -K-HfiT) + 7r-iAarg(L(s, /)). 

Since L{s,f) 7^ if ^(s) — 4 > 1/2, Aarg(L(s,/)) equals the variation of the 
argument of L{s, /) along the line segments Li = [4, 4 + i + ^u] , L2 = [4 + i + ^u, 4 + 
l+fi + iT], L3 = [4 + ^ + fi + iT,4 + iT] for any ^u > 0. Take T = 100, ^ = 1. Then 
we have 7r~^'i?/(100) = 69.0171 ■ ■ ■ . Hence if we can show |Aarg(L(s, /))| < tt/2, we 
can conclude N(T) = 69. For this purpose, it suffices to show that K(L(s, /)) > 
when s G Lj, z = 1, 2, 3. For Li, we have shown this fact above. For L2, this 
fact can be proved as in [Y2], §4. For L3, we divide it into 150 small intervals and 
appeal to our heuristic calculation. We have observed 

3f?(L(s,/))>0.83 on L3. 

Thus we conclude that the Riemann hypothesis holds for L(s, /) when < Q{s) < 
100. All the zeros are simple. 

Now let us consider zeros of L(s, g). We have g{iy) > for y > since 9{iy) > 0, 
r]{iy) > for 2/ > 0. By the integral representation 



{2n)-^T{s)L{s,g)= / g{iy)y^-Uy, 
Jo 



we see that L(s, g) > for s > 0. For T > 0, let N{T) denote the number of zeros 
of L(s, (7) counted with multiplicity in the domain^ 

|3f?(s)-9/4| <2, 0<S(s)<T. 

By (1.7) taking 5 = 2, we have 

N{T) = TT-'MT) + 7r-iAarg(L(s, g)). 

Here Aarg(L(s, (7)) denotes the variation of the argument of L{s,g) along the line 
segments Li = [9/4,17/4], L2 = [17/4, 17/4 + iT], L3 = [17/4 + zT, 9/4 + zT]. 
Take T = 100. Then we have 7r-i't?g(100) = 79.1885 ■ ■ ■ . Dividing L2 and L3 into 
small intervals, we have observed 7r~^A(argL(s, (/)) = 0.8114 ■■■. Thus we obtain 
N{T) = 80. On the otherhand, we have obtained only 54 zeros on the critical line. 
Therefore, assuming that these zeros are simple, there must exist 13 zeros in the 
right-hand side of the critical line: 9/4 < ?R:{s) < 17/4, < S(s) < 100. These 
zeros, together with those in 100 < Q'(s) < 150, are given in Table 3.5. 

Our method of calculation of these exceptional zeros is as follows. Let us consider 
a box B given by di < ^(s) < d2, hi < $5(s) < /i2. By the argument principle, we 



^The zero free region of L{s,g) is non-trivial. Here we content ourselves by regarding \^{s) 
9/4| < 2 is "sufficiently wide". 



can determine whether L{s,g) has a zero incide B or not. First we find a box B 
in which L(s, g) has zeros by trial and error. Then dividing B into sub-boxes and 
applying the principle above successively, we can obtain a good approximation for 
a zero inside of B.'^ 

The Riemann hypothesis does not hold for L(s, g). This should be of no surprise 
since L{s,g) does not have an Euler product. 

We shall study one more example of modular forms of half integral weight. For 
/c > 1, put 

n=l d\n,d>0 

Let 

^(z) = 

2Txi 



fin 

5{z) = {2G^{Az)e'{z) - G'^{Az)e{z)) = Yl ^(^)^''- 



n=l,n=0,l mod 4 



Then we have b G 5'i3/2(ro(4)) (cf. Kohnen-Zagier [KZ], p. 177) and b corresponds 
to A(2) = ?7(z)^'^ G S'i2(S'L(2, Z)) under the Shimura correspondence. The values 
of c(n) can easily be computed by 



[V] 



c(n) = u){\fn) ■ n + 120 >^ uj{\/n — Am)a^{rn){2n — 9m) — lbna^{n/A), 
where 

r 1 if X G z, 

^ ' 1 if X ^ Z. 

Let 

5o{z) = Yl cin)q-/\ 

n=0 mod 4 

Then we have do G S'i3/2(ro(4)) and 

Ay 
(cf. [KZ], p. 190.) Put 

h±{z) = (5t2-^(5o)(^), R{s,h±) = 2^(27r)-T(s)L(s,/i±). 

Then we have 

/i±(f) = ±(2y)6+i/2/i±(zy), y>0. 

Hence we obtain the functional equations for the entire functions R{s, h±): 

R{—-s,h±) = ±R{s,h±). 



^It is more efficient to use a variant of the Newton method once we get a rough approximation. 
After finding a precise location, final check should be done by the method described above. 
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We have computed zeros of L(s, h±) on the critical hne 3?(s) = 13/4 in the range 
< Q'(s) < 100. The results are given in Table 3.6. There u^ denotes the n-th 
zero of L(s, h±) for s = ^ + iu. 

In Table 3.7, we also list zeros of L{s,h±) not on the critical line which are 
searched in the region 13/4 < K(s) < 25/4, < $5(s) < 100; p^ denotes a zero of 
L(s,h±). It is a very interesting phenomenon that L(s,h-) has much more zeros 
outside of the critical line compared with L(s, h+). 

§4. L-functions attached to Hecke characters of infinite order of real quadratic fields 

In this and the next section, we shall study two types of L-functions which 
are closely related to algebraic number fields. We can still apply our method of 
calculation described in §1 efficiently. However the situation changes drastically. 
The repeated application of partial summation does not yield good results beyond 
rather limited number of times. Thus our calculation cannot be as accurate as in 
the case of modular forms treated in §2 and §3. 

Let k he a, real quadratic field. For simplicity, we assume that the class number 
of /c is 1. Let D be the discriminant and e be the fundamental unit of k. Let k^ 
denotes the idele group of k. For a finite place v of /c, let ky denote the completion 
oi k aX V and Oy denote the ring of integers of ky. Since k is of class number 1, we 
have 

(4.1) fc^ = A;^(JjD^ xR^ xR^) 

V 

where v extends over all finite places of k. Let x = Ylv Xv be an unramified unitary 
Hecke character of k^. Let ai (resp. (72) be the identical (resp. non-identical) 
isomorphism of k into R and let ooi (resp. 002) be the corresponding archimedean 
place of k. As unitary characters of R^, Xooi and X002 take the following form: 

(4.2) XooA^)=^Sn{xr^xr for a:eK^^. -R^ j = l,2, 

where rrij = or 1, Vj G R. By (4.1), we see that x is completely determined 
by Xooi and Xoo2- Since x is trivial on /c^, we must have x(x) = 1 for all k^ fl 
(Yly Oy X R^ X R^), which is the group of units of k. Therefore we have 

(4.3) (_i)"^i+"^2 = i^ sgn(e'^i)"^isgn(e'^2)"^2|e'^i|*^i|e'^2p"2 = 1. 

It is easy to see that (4.3) is a necessary and sufficient condition for x, which is 
determined by Xooi and X0027 to be a Hecke character of k^. By (4.3), we have 
mi = 1712- Put m — mi. Then (4.3) is equivalent to 

(4.4) |e|*(«i-«2) = sgn(iV(e))"". 

Let X* be the associated ideal character of k. If (a), a E k^ is a prime ideal, we 
have, by definition 



x((«-' 



X* 


^((«)) = 


-x{{- 




,1,-- 


■ ,a, ■ 


■■,1 


J 


■)) 










,1, 




■ ,a" 


-1 . . , 

7 ' ' ' 


■)) 


= sgn(A^(a))"" ■ 


i\a^^ 


ivi 


a"^ 


\iv2\- 


-1 
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Here (■ ■ ■ , 1, ■ ■ ■ , a, ■ ■ ■ , 1, ■ ■ ■ ) G /c^ denotes the idele whose (a)-component is a 
and aU the other components are 1. Hence we have 

(4.5) x*((a)) = sgn(A^(a))""|a'"ir'^^|a'^2r™2 for every a e A;^, 
L(s, x) = L{s, X*) = E(«) X*((a))iV((a))-^ Put 

R{s,x) = |I)|'/V-("+"^)7r-*("i+"2)/2r((s + m + zt;i)/2) 

r{{s + m + iv2)/2)L{s,x), 
R{s,x~^) = |I5|'/'7r-(^+"^V*("i+"2)/2r((s + m-it;i)/2) 

T{{s + m-iv2)/2)L{s,x~^)- 
Then the functional equation is (cf. Weil [W], Langlands [LL]) 

(4.6) R{s, x) = {-irx*{{d))R{l - s, X-'), 



where (d) denotes the different of k over Q. Since R{s, x) = R{s, X ^)5 we can put 
(4.6) in the form of (1.2): 



(4.7) R{1 -s,x) = {-irx*{{d))R{s, x)- 

We get 

(4.8) x*m)-'/^Ris,x) = {-irx*mr'^^Ris,x) if 3?(s) = i/2. 

Hence x*(((i))~^'^i?(s, x) takes real or pure imaginary values on the critical line 
according as to = or 1. We also note that we may assume t'2 = without losing 
any generality since the choice of V2 can be taken into account as the shift of the 
variable s. Then, if to = 0, we have vi = — 2?i7r/loge with n G Z by (4.4). We 
denote this Hecke character by Xn- By (4.5), we have 

(4.9) (;^^),((a)) = |a|2--/iog.^ « ^ fc^ 

If m = 1, by (4.4), we have vi = —2mi/\oge or vi = —{2n + l)7r/loge with n G Z 
according as A^(e) = 1 or A^(e) = —1. We denote this Hecke character by Xn- By 
(4.5), we have 

.4... r sgn(iV(a))|«|2Wioge if iV(e) = l, 

^ ' ^^^>*^^''>> I sgn(iV(a))|a|(2-+iWiog. jf N{e) = -1. 

As our first example, we take k = Q(a/2). We have e = y/2 + 1, (d) = (a/2)"^, 
\D\ = 8. We are going to study L(s, xi) applying our summation method. For 
s = I + zt, t = 15 and 50, the values of 

Rj = ^{x*mr'/^expiiMt))si^^). ij = '^{x.mr''^^MiMt))s^N) 

are given in Tables 4.1a and in 4.1b respectively. Here X*(((i))~^'^ = v2 , 

^o{t) = arg(8^/27r-%27ri/iog£p^g _ 27rz/loge)/2)r(s/2)). From this table, it is 
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evident that R/s for higher j do not give good results. We can judge, from the 
values of |/j|, i?2 gives the best result, then i?3, Ri in this order. We must be 
more cautious than in §2 and §3 about the accuracy of the value e*'^°*^*^L(s, xi). 
For example, let t = 15. We empirically judge that e*^"*^*-'L(s, xi) = 2.17375 with 
error ^ 10~^ from R2 and l2- We have constructed Table 4.2 in which zeros on 
the critical line K(s) = 1/2 are listed in the range \'^{s)\ < 30. Here, for n > 1, t^ 
(resp. t-n) denotes the n-th zero of L(s,xi), s = 1/2 + it on the critical line for 
< t < 30 (resp. < -t < 30). 

Let us examine the Riemann hypothesis for L(s, xi) in the range < '^{s) < 100. 
We have observed 84 sign changes of e''^°(*^L(i +zt, Xi) for < t < 100. For T > 0, 
let N{T) denote the number of zeros of L(s, xi) counted with multiplicity in the 
domain 

< ^{s) < 1, < 53(s) < T. 

Taking S = 1/2 in (1.7) , we get 

(4.11) N{T) = 7v-^^{T) + 7r-iAarg(L(s, xi))- 

Since L(s, xi) 7^ if ^(s) > 1, Aarg(L(s, xi)) equals the variation of the argument 
of L(s,xi) along the line segments Li = [i, 1 + (u], L2 = [1 + n,l + fi + iT], 
L3 = [1 + fi + iT,l + iT] for any ^u > 0. Take T = 100 and fi = 1. We have 
7r-^t?(100) = 84.8864 ■••. Hence if we can show -n < Aarg(L(s, /)) < 0, then 
we can conclude that A^(IOO) = 84. For Li, we divide it into 15 intervals. We 
observed that L{s,xi) moves from 0.3482 + 0.0712v/^ to 0.8011 + 0m69^/^ 
keeping 3f?(L(s, xi)) > 0. For L2, we can show without difficulty that 3f?(L(s, xi)) > 
on L2. For L3, we divide it into 150 small intervals. We observed that L(s,xi) 
moves from 0.8159 + 0.1227v/^ to -0.0110 - 0.0072 v^^ when s moves from 2 + 
100^/^ to i + lOOv/^; L(s,xi) never crossed the half line Q'(L(s, xi)) = 0, 
K(L(s, xi)) ^ 0. Hence we get A^(IOO) = 84. The Riemann hypothesis holds and 
all zeros of L(s, xi) are simple zeros in this range. 
By (4.11), we should have 

712 72 
A^(IOO) = 84.8864 tt"^ arctan( ) - (1 - n'^ arctan( )) = 84.0067 ■ ■ ■ . 

The error is about 6.7 x 10~^ and this is much bigger than the usual error inherent 
in our calculations. The reason is that L(s,xi) takes rather small value at | + 
100a/— 1 ; such an error can be made much smaller in the following way. We take 
T — 101. We have observed 86 sign changes of e'^'^°^^^L{^ + it, xi) for < t < 101. 
We have 7r-i^?(101) = 86.0881 ■ ■ ■ . We divide [2 + lOlv^, ^ + lOlv^^] into 150 
small intervals. We observed that L(s, xi) rnoves from 0.9322 + 0.2563V— 1 to 
1.7137 - 0.1288 V^. Hence by (4.11), we have 

71 2 1 288 

iV(lOl) = 86.0881 TT"^ arctan( ) - tt"^ arctan( ) = 86 

^ ^ ^3482^ 47137^ 

with error less than 10""^. 

In Table 4.3, we have listed 44 examples of L(s, Xn) for which we made experi- 
ments in the range < Q'(s) < T; N{T) denotes the number of zeros of L{s, Xn) in 
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the domain < K(s) < 1, < Q'(s) < T. We found that aU zeros in the ranges of 
Table 4.3 he on the critical line and are simple. 

§5. Artin L-functions 

Let k be the minimal splitting field of the irreducible polynomial f{X) = X^ — 
X+1 over Q. Then k D Q(Vl9 ■ 151) = /cq, k is unramified over /cq, Gal(/c/Q) = 5*5, 
Gal{k/ko) = A^. The discriminant A of a root of f{X) is 19 ■ 151. This example 
is due to E. Artin (cf. Lang [LG], p. 121). Let p be an irreducible 4-dimensional 
representation of 5*5 whose character Xp is given as follows. 

conjugacy class (1) (12) (123) (1234) (12)(34) (12)(345) (12345) 
Xp 4 2 10 -1-1 

Since 5*5 does not have a subgroup of index 4, p is not monomial. We have L(s, p) = 
Yip Lp{s, p) for 3f?(s) > 1 with the Euler p-factor Lp(s, p). We can compute Lp{s, p) 
as follows. First we assume that a prime number p is unramified in k, i.e., p ^ 19, 
151. Then we see easily that Lp(s, p)~^ equals 

(l_p-)3(i+p-.) if ap={(l2)}, 
(l_p-)2(i+p-+p-2s) if ap={(123)}, 
(l_p-)(l+p-)(l+p-2^) if a^={(l234)}, 

l+p-'+p-'^'+p-'^' +p-^' if (7p={(12345)}, 
(l_p-)2(i+p-)2 if ^^={(i2)(34)}, 

(l-p-^)(l+p-^)(l+p-^+p-2^) if (7p={(12)(345)}. 

Here Up denotes the Frobenius conjugacy class of p and {r} denotes the conjugacy 
class of r G 5*5. 

Let p = 19 or 151. Let /p denote the inertia group of a prime factor p of p in k. 
By definition, we have 

(5.1) Lp{s,p)-^ = det(l - {p{a,)\V''')-p-^). 

Here V denotes the representation space of p, F^p the subspace of /p -fixed vectors 
and (Tp a Frobenius of p which is determined modulo /p. Since k is unramified 
over /co, it is obvious that \Ip\ — 2, Ip <^ Gal(/c//co) — ^s- Hence we may assume 
that /p is generated by (12) choosing a suitable p lying over p. Let Dp denote the 
decomposition group of p. Then Dp > /p and Dp/ Ip is generated by ap mod /p. 
We have 

iVs3(/p)=/pX53, 

where A^5g(/p) denotes the normalizer of /p in 5*5 and 5*3 denotes the permutation 
group on three letters {3, 4, 5}. 
Let p = 19. Then 

f{X) = {X -6)^{X^ + 12X'^ + 13X + 9) mod 19 
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is the factorization of f{X) mod p into irreducible factors in (Z/pZ) [X] . Therefore 
the residue field extension Dfc/p of Z/pZ contains the cubic extension of ZjpZ, 
where Dfc denotes the ring of integers of k. Hence we immediately obtain D^/Ip = 
Z/3Z and that cxp may be taken as (345) G 5*3 . Now we find easily that 

Lp(s,p)-^ = l-p-3^ if p = 19. 

Let p ^ 151. Then 

f{X) = {X- 39)2(X - 9)(X2 + 87X + 61) mod 151 

is the factorization of f{X) mod p into irreducible factors in {Z/pZ)[X]. By a 
similar consideration as above, we find that Dp/Ip = Z/2Z and that ap may be 
taken as (34) e 5*3 . We obtain 

L^{s,p)-^ = {l+p-^){l-p-'f if p = 151. 

Let /(p) denote the Artin conductor of p. We easily obtain 

/(p) = 19-151. 



rsj 



For example, let p = 19 and p be as above. We have shown Gal(/cp/Qi9) = Dp 
Z/2Z X Z/3Z. We find that the restriction of p to -Dp splits into a direct sum of 
four one dimensional representations of Dp such that three of them are unramified 
and one is ramified. Hence the exponent of 19 in f{p) is 1. 

We take an isomorphism a of /c into C and let c G Gal(/c/Q) = 5*5 be the 
restriction of the complex conjugation to k. Then c G Gal(A;//co) — ^5- Hence c is 
conjugate to (12)(34) in 5*5. Let Gal(C/R) be identified with the decomposition 
group (c) of the archimedean place of k which corresponds to a. The restriction 
of p to Gal(C/R) splits into a direct sum of two trivial representations and two 
non-trivial representations. Therefore the Gamma factor to go with L{s, p) is given 
by (cf. Langlands [LL]) 

(7r-^/2r(s/2))2(7r-(^+i)/2r((s + l)/2))2. 

Put 

R{s, p) = (19 ■ l^lY/\-^'T{s/2fT{{s + l)/2fL{s, p). 

Since p is equivalent to its contragredient, we have the functional equation 

(5.2) R{s,p) = tiR{l-s,p), 

where k = ±1 is the Artin root number attached to p. Let if) be the additive 
character of Q^/Q such that 

^l^ooix) = exp(27rV^a:), x G Qoo = R, 

ipp{x) = exp(— 27rv^^Fr(x)), x G Qp, 

where Fr denotes the fractional part of x. By a theorem of Langlands, we have 

V 
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with the e- factor defined in [LL]. By the above considerations, we easily get 



1 / xf 1 if P 7^ 19, 151, 

2'^^'"^^^ \i if p = 19 or 151, 



Hence we obtain 

(5.3) K^l. 

The values of 

Rj = 3?(exp(zz?(t))>5}^^), Ij = S(exp(zz?(t))5j^^) 

for L(s, p), s = I + zt, t = 5 are given in Table 5.1. Here 

i}{t) = arg((19 ■ 15iy/\-^'T{s/2fT{{s + l)/2)2), g = - + it. 

In Table 5.2, we give the values of -u^, the n-th zero of L{s^p), s = | + m on the 
critical line for < tt < 10. 

§6. Estimation of errors in our calculations 

The most serious defect of our method of calculation is that we do not have rig- 
orous controle of error estimates. In previous sections, we regarded the magnitude 
of '^{e^'^^^' S j^) (resp. 'R{e'^'^^^' S j^) ) as a rough measure of errors from the true 
value, when e^"^^^' L{a + it) should be real (resp. pure imaginary). In this section, 
we shall present several data which support this practice. 

Suppose that the functional equation (1.2) for L{s) holds. Then we have 
e^'?(*)L(/c/2 + it) gR, t G R where z?(t) = arg(KiA^^ J]^^ r(6iS + q)), s = k/2 + it 
in the notation of (1.4). Take < ti < t2 so that e*''*^*-' rotates on the unit circle 
exactly once when t moves from ti to t2- We expect that maxt^<t<t2 {^{e'^'^^^^Sj^ )\ 
can be used as the measure of errors. More explicitly, it seems plausible that 

(6.1) \^{e'^^'^S^^^)-L{k/2 + it)\<10 max |3(e^^(*)5{J^)|. 

ti ^t^t2 

In examples below, we use ^{S]^ ) as a substitute for L{k/2 + it) taking large M 
and p (except for in Example 6), and examine the ratio of two terms in (6.1) for 
Sj^ taking relatively small A^ and /. The results are given in Table 6.1. 

Example 1. We take the primitive form / G 5*8 (Fq (2)) and consider the L- 
function L(s, /), s = 4 + it. When t moves from 97.9 to 100, e *^*-' rotates on the 
unit circle approximately once. We calculated the ratio 

1) _ max97.9<t<ioo(|»t(e ^^^20oo~^ '^looooJIJ 



^5Af 



max97.9<t<ioo(|3(e'''(*)5'^ooo)l) 
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for 1 < A^ < 6, dividing [97.9, 100] into 21 intervals of length 0.1. 

Example 2. We take g G 5*9/2 (Fq (4)) as in §3. When t moves from 98.2 to 100, 
^i-'&i't) rotates on the unit circle approximately once. We calculated the ratio 

^(2) ^ max98.2<t<ioo(|g^(e--^(^)4oi'J ' e--^^^^>gSo)l) 
max98.2<t<ioo(|S(e^'*(*^5'2S)|) 

for 1 < A^ < 6, dividing [98.2, 100] into 18 intervals of length 0.1. 

Example 3. We take A G 5'i2(5'L(2, Z)) and consider the L-function L(^)(s, A), 
s = 17 + it, attached to the third symmetric power representation of GL{2). When 
t moves from 17.5 to 20, e"^"^^^' rotates on the unit circle approximately once. We 
calculated the ratio 

j3) _ maxi7.5<t<2oU-^ie ^^52000"^ '^loooi 



„. , _ - - .. . ^^^„ aoooo^ 

' 5N ~ 



niaxi7.5<t<2o 



ime^^wsg^^)\) 



for 1 < A^ < 6, dividing [17.5,20] into 25 intervals of length 0.1. In this example, 
we have normalized z?(t) as in §2 so that e *^*''L*^"^''(17 + it, A) is pure imaginary. 

Example 4. We consider the L-function L^'^\s,A), s = 45/2 + it, attached to 
the fourth symmetric power representation of GL{2). When t moves from 7.2 to 
10, e^"^^^' rotates on the unit circle approximately once. We calculated the ratio 

^(4) ^ maxr.2<t<io(|g^(e--^(^)4^o^J - e^'^i'^ S[lll,)\) 
'"^ maxT.2<t<io(|S(e'''W4S)|) 

for 1 < A^ < 6, dividing [7.2, 10] into 28 intervals of length 0.1. 

Example 5. We consider the Artin L-function treated in §5. When t moves from 
4 to 5.8, e*'*'^*-* rotates on the unit circle approximately once. We calculated the 
ratio 

^(5) ^ max4<t<5.8(|g^(e--^W>ggSo - e^'^'^ Sg^jjl) 

max4<,<5.8(|S(e^''W5gSo)l) 
for < A^ < 5, dividing [4, 5.8] into 18 intervals of length 0.1. 

Example 6. We consider the Hecke L-function L{s, xi) for k = Q(v^) treated in 
§4. When t moves from 47.9 to 50, e*^'^*-* rotates on the unit circle approximately 
once. We calculated the ratio 

^(6) ^ max47.9<t<5o(|g?(e--^(*)giSoo - e^-^^^^^gLo)!) 
max47.9<t<5o ( I $5(e*'^(*) slgfj^g) \ ) 

for < A^ < 5, dividing [47.9, 50] into 21 intervals of length 0.1. 

§7. A comparison with the explicit formula 

We take the new form / G 5'8(ro(2)) treated in §3. Let n be the irreducible 
unitary automorphic representation of GL{2, Qa) which corresponds to /. We 

have 

7 
Lf{s,w) = L{s + -,f), 
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where Lf{sj tt) denotes the finite part of the Jacquet-Langlands L-function attached 

to TT. Let 

p 

be the Euler product of Lf{s, tt). For p = 2, the Euler 2-factor degenerates so that 
a2 = — 1/v^, (32 = 0. For p 7^ 2, we have \ap\ = \(3p\ = 1 by the Ramanujan- 
Petersson conjecture proved by P. Dehgne. We have the exphcit formula 



E E K+/3;)logP = - lim Yl ^ ^' 



L'Mtt) 



(7 2) P l<n,p"<x \Q(p)\<T ' /v 5 y 

+ log(^^) - 2X-V2 _ 2-3/2 _ 2 _5/2 
V2; — 1 3 5 

for a; > 1. Here J2 means that the term (a^ + (3p) logp should be multiplied by 
1/2 when p"^ = x; p extends over zeros of Lf{s, it) such that < 3f?(p) < 1. This 
formula can be shown in the usual way as in Ingham [In], p. 77-80. The last term 

9{x) := log{y^) - 2x- V2 _ 2-3/2 _ 2-5/2^ 

which is equal to 

^E 2A; + 1 ' 

fc=3 

represents the contribution of the trivial zeros of Lf{s, tt); they are at s = — |, — |, 
_ii 

2 ' ■ ■ ■ • 

Now it seems very interesting to compare both sides of (7.1) numerically using 

the zeros of L(s, /) given in Table 3.3. We approximate 

lim 7 — 

T^+oo_ ^ _ p 

by 

(7.2) h69{x) := Y^ \fx{t\ + -)-^ [cos(t„ log x) + 2t„ sin(t„ log 



mP)\<T ^ 



with tn given in Table 3.3. We have 

7 7 

L/(0, tt) = L(-, /) = 0.5942254156 ■ ■ ■ , L'.(0, tt) = L'(-, /) = 0.1875716234 ■ ■ ■ , 

2 ■'2 

L'(0,7r) 

-^^ = 0.3156573558 ■ ■ ■ . 

L/(0,7r) 

To obtain these values, we simply applied our repeated abel summation technique 
as before, though more rigorous evaluation could be made in this case. 
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In Figure 7.1, we drawed the graphs of the "step function" 

E E' K+/3;)i°gp 

P l<n,p"<x 

and —hQg{x) + g{x) — 0.3156 for 1.1 < x < 20. (We have used "Mathematica" to 
make Figure 7.1.) The coincidence seems fine. 

§8. Sample programs 

In this section, we shaU present a few sample programs to compute L(s, /), 
/ G 5'8(ro(2)) (cf. §3). All programs, which are ready to be executed, are written 
in UBASIC created by Y. Kida. Let 

oo 

f{z) = {v{zM2z)f = J2^r.qr 

n=l 

Using Program A, we can compute a^ for 1 < n < M for any M, 1 < M < 
lO'*. From line 50 to 150, the coefficients A{n) in ri{z) = Q'^ Xl^i A{n)q^~^ are 
computed for 1 < n < M using Euler's formula 



v{z)=q'/''ll(l-qn = Q'^'' E ("1)'^^^ 



n(3n+l)/2 

-1- ; k' 

n=l n= — oo 



From line 160 to 220, the coefficients B(n) in rj{z)r){2z) = q^'^ Xl^i B{n)q^ ^ are 
computed for 1 < n < M by B{I) = Ej+2L=7+2 MJ)ML)- From hue 270 to 390, 
the expansion of r]{z)ri{2z) is raised to the eighth power; the final result will be 
stored in the data file "wt8 " . 

Program B computes the value of e*'^^'^*^L(4 + zt, /) for t = 100. To save the 

space, this program gives the values of e'^^^'''^'^SjJ for A^ = 2000, < / < 40. By 
point 15 command (line 20), UBASIC gives the precision to the 70-th digit. From 
line 300 to 580, the value of e^'^f^^' will be computed and stored in the variable 
THE. The calculation proceeds as follows. We have 

^/(t) = arg(2^/2(27r)-T(s)),=4+^t = {\ log2 - log27r)t + argr(4 + it). 
We also have 
argr(2) = argr(2 + 1) - arctan($5(2)/3f?(z)), ^{z) > 0, argr(z) = 53(logr(z)). 

By these formulas, it suffices to compute log r{z + 100) for ^ = 4 + it. We have the 
asymptotic expansion (cf. [WW], p. 252) 

1,^ 1 >^ {-ly-^Br 



U) logT{z) r^ {z - -)logz - z + -log2n + J2 



-^ 2r(2r-l)z2'^-i 



where B^ denotes the r-th Bernoulli number. In (8.1), we use the terms up to 
r = 10. We can get an approximation of e f^^' which is accurate at least to the 
40-th digit (cf. [WW], p. 252). 
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From line 660 to 740, the values of 



N N 



4°) = 5:a.n- .« = $:4'-^), l</<40, r«('^) 

n=l n=l 



^n ^ri) 



are computed. From line 770 to 1020, the values of Sp^ are computed for 1 < / < 40 
using (1.10) and (1.11). We use the approximation 

100 I 



(0 






The error from the truncation by 100 is negligible (cf. §1). The variable U{K) 
stands for (_i)fc£(£±ihi£±^^. the variable Co stands for E!n=i(-l)"'(ri)^^ 
the variable T{I) stands for 



100 I 



I — ; m — 1 \ / 



A:=Z m=l 



In the line 1000, X = S^S^ is multiplied by e''^i^^\ 
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Program A 



10 word 8 

20 point 2 

30 dim A(10000),B(10000) 

40 input M 

50 Al=sqrt(24*M+l) 

60 Ml=int((Al+l)/6) 

70 A(l)=l 

80 for 1=1 to Ml 

90 Il=I-2>Kint(I/2) 
100 12=1-2*11 
110 J=int((3*M+I)/2)+l 
120 A(J)=A(J)+I2 
130 J=int((3*M-I)/2)+l 
140 A(J)=A(J)+I2 
150 next I 
160 for 1=1 to M 
170 Il=int((I+l)/2) 
180 for L=l to II 
190 J=I-2>kL+2 
200 B(I)=B(I)+A(J)>kA(L) 
210 next L 
220 next I 
230 for 1=1 to M 
240 A(I)=B(I) 
250 B(I)=0 
260 next I 
270 K=l 
280 for 1=1 to M 
290 print K,I 
300 for J=l to I 
310 B(I)=B(I)+A(J)*A(I+1-J) 
320 next J 
330 next I 
340 for 1=1 to M 
350 A(I)=B(I) 
360 B(I)=0 
370 next I 
380 K=K+1 

390 if K<4 then goto 280 
400 open "wt8" for output as [Jl 
410 for 1=1 to M 
420 print I,A(I) 
430 print 91,A(I) 
440 next I 
450 close p 
460 end 
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Program B 

10 word 70 

20 point 15 

30 dim Bn(20),Bd(20),C(2000),T(100),U(200),Sm(50) 

40 Ab=40 



50 


Bn(l): 


=1 


60 


Bd(l> 


=6 


70 


Bn(2> 


=1 


80 


Bd(2> 


=30 


90 


Bn(3> 


=1 


100 


Bd(3) = 


=42 


110 


Bn(4) = 


=1 


120 


Bd(4) = 


=30 


130 


Bn(5> 


=5 


140 


Bd(5) = 


=66 


150 


Bn(6) = 


=691 


160 


Bd(6> 


=2730 


170 


Bn(7> 


=7 


180 


Bd(7) = 


=6 


190 


Bn(8) = 


=3617 


200 


Bd(8> 


=510 


210 


Bn(9> 


=43867 


220 


Bd(9) = 


=798 


230 


Bn(10 


)=174611 


240 


Bd(10 


)=330 



250 A=l/sqrt(3) 

260 P=6*atan(A) 

270 S=4+100*rti 

280 Ss=S 

290 T=im(S) 

300 Th=T>K(log(2)/2-log(2>KP)) 

310 U=0 

320 for 1=1 to 100 

330 Sl=re(S) 

340 S2=im(S) 

350 if S1>S2 goto 380 

360 U=U-(P/2)+atan(Sl/S2) 

370 goto 390 

380 U=U-atan(S2/Sl) 

390 S=S+1 

400 next I 

410 Th=Th+U 

420 T=im(S) 

430 R=re(S) 

440 Xl=atan(T/R) 

450 X2=(log(R*R+T*T))/2 

460 X3=X2+Xl*Hi 

470 X4=(S-(1/2))*X3 
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480 Th=Th+im(X4)-T 

490 K=10 

500 S1=S 

510 for 1=1 to K 

520 Xl=Bn(I)/(2*I*(2*I-l)*Bd(I)*Sl) 

530 Th=Th+im(Xl) 

540 S1=-S1*S*S 

550 next I 

560 Xl=int(Th/(2*P)) 

570 X2=Th-2*P*Xl 

580 The=exp(X2*tti) 

590 S=Ss 

600 X=0 

610 open "wt8" for input as jjl 

620 for M=l to 2000 

630 input tll,C(M) 

640 next M 

650 close IJl 

660 for M=l to 2000 

670 Sm(l)=Sm(l)+C(M) 

680 for 1=2 to Ab 

690 Sm(I)=Sm(I)+Sm(I-l) 

700 next I 

710 Xl=log(M) 

720 X2=exp(-Xl*S) 

730 X=X+C(M)*X2 

740 next M 

750 Z=The*X 

760 print Z 

770 N=2001 

780 A=log(N) 

790 Nl=exp(-A*S) 

800 U(l)=-S 

810 for K=2 to 100 

820 U(K)=-U(K-1)*(S+K-1)/K 

830 next K 

840 for 1=1 to Ab 

850 T(I)=0 

860 if I>1 then goto 890 

870 X=X-Sm(l)*Nl 

880 goto 1000 

890 Ii=I-l 

900 for K=Ii to 100 

910 Co=-Ii 

920 Al=-Ii 

930 for L=2 to li 

940 Al=-Al*(Ii-L+l)/L 

950 Co=Co+(Al*(LAK)) 

960 next L 
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970 T(I)=T(I)+Co*U(K)/(NAK) 

980 next K 

990 X=X-Sm(I)*Nl*T(I) 

1000 Z=The*X 

1010 print I,Ss,Z 

1020 next I 

1030 end 



^9. Conjectures 



In this section, we shall discuss a few conjectures which emerged during the 
process of experiments: No non-trivial coincidences of zeros of two L-functions 
attached to two non-equivalent irreducible A-adic representations of Gal(Q/Q) are 
found. We shall use the framework of automorphic representations of GL{n^ Qa) 
to formulate this fact in more general case. 

Let TT be an irreducible unitary cuspidal automorphic representation of 
GL{n, Qa). The contragredient representation tt to tt is equivalent to the complex 
conjugate representation yf of tt and we have the functional equations: 



L(s, tt) = e(s, 7r)L(l — s, tt), L(s, tt) = L(s, tt). 

Let WjT be the central character of tt. For s G C, set f^i^x) = |a;|^, x G Q^ where 
\x\a denotes the idele norm of x. We can find a t G R so that Wtt'^'* is a character 
of Q^ of finite order. Since Uj^u'^^ is the central character of tt ® (//**/"■ o det) and 
L{s + —,Tx) = L(s, TT (g) (z/**/" o det)), we may assume, without losing substantial 
generality, that cUjr is of finite order. 

Conjecture 9.1. Let tti and 1x2 he irreducible unitary cuspidal automorphic repre- 
sentations of GL{n,(^A) with the central characters Ut^^ andui-,^^. We assume that 
TTi is not equivalent to 112 and that 0;,^ and cu^ra ctre of finite order. Then L(s, tti) 
and L{s, 712) have no common zeros in the critical strip < K(s) < 1 except for 
s = 1/2. 

Remark. If we replace Q by an algebraic number field, the assertion is obviously 
false. 

As a variant of 9.1, we can formulate a conjecture on L-functions of motives. Let 
E be an algebraic number field of finite degree. Let Mi and M2 be motives over Q 
with coefficients in E of pure weights t^i and W2, of ranks ni and n2 respectively. 
We assume wi = W2 and put w = wi. Fix an embedding a of E into C and let 
L(s,Mi) (resp. L(s,M2)) be the L-function of Mi (resp. M2) with respect to a. 
For a finite place A of E, let pi : Gal(Q/Q) — > GL{ni, Ex), i = 1, 2 be the A-adic 
representation obtained from the A-adic realization of Mj. We assume meromorphic 
continuation of L(s, Mi), i = 1, 2 to the whole complex plane. 

Conjecture 9.2. Assume that pi, i = 1, 2 are absolutely irreducible and that pi 
is not equivalent to p2 for a finite place A of E. Then L(s, Mi) and L(s, M2) have 
no common zeros in the critical strip ^ < ^(s) < ^ + 1 except for s = {w + l)/2. 

Remark. We understand that a pole of order k is a, zero of order —k. Conjecture 
9.2 implies the (usual) Artin conjecture except for s — 1/2. 
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Remark. The condition on pi implies that Mi and M2 are simple motives with 
coefficients in Q. D. Blasius observed that the analogous conjecture for motives over 
a finite field is true under the Tate conjecture (cf. Milne [M], p. 415, Proposition 
2.6). 
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Table 2.1 






A^ 


-^0 -Ro 


h 


-R5 iio 


-Rio 


1000 


3.95 -9.3 X 10"^ 


^ 3.9415443 


1.3 X 10-3 3.9237277 


3.0 X 


10-3 


2000 


4.02 3.3 X 10"^ 


^ 3.9248413 


4.6 X 10-3 3.9238190 


-1.1 X 


10-5 


4000 


3.90 1.0 X 10"- 


^ 3.9233635 


-2.1 X 10-4 3.9242209 


-1.8 X 


10-5 


6000 


3.88 -3.9 X 10"^ 


5 3.9239218 


-3.6 X 10-4 3.9242178 


-3.7 X 


10-6 


8000 


3.85 -1.4 X 10"^ 


^ 3.9245956 


6.6 X 10-5 3.9241989 


5.7 X 


10-6 


10000 


3.92 2.5 X 10"^ 


^ 3.9240769 


-8.5 X 10-^ 3.9242054 


-3.2 X 


10-^ 














A^ 


/l5 it 


'15 


ho R20 




1000 


3.92361007 -1.3 


X 10-4 3.924178604 -4.1 x 10"^ 




2000 


3.92423322 -3.8 


X 10-5 3.924207205 4.9 x 10"^ 






4000 


3.92420427 1.8 


X 10-6 3.924203509 -1.8 x 10-^ 






6000 


3.92420338 4.6 


X 10-^ 3.924203748 -2.9 x 10-^ 






8000 


3.92420375 -1.9 


X 10-^ 3.924203739 7.2 x 10-^ 






10000 


3.92420370 2.4 


X 10-s 3.924203738 -1.0 x 10-^ 
















N 


/25 


-R25 


ho R30 






1000 


3.92438069843 


8.8 X 10-6 


3.92417001558 2.0 x 10" 


-4 




2000 


3.92419997465 


1.2 X 10-6 


3.92420393248 -9.9 x 10" 


-7 




4000 


3.92420378031 


2.6 X 10-s 


3.92420373301 3.0 x 10" 


-9 




6000 


3.92420373864 


1.5 X 10-9 


3.92420373759 -1.7 x 10" 


-10 




8000 


3.92420373800 - 


-2.7 X 10-1° 


3.92420373818 -8.9 x 10" 


-13 




10000 


3.92420373814 


5.8 X 10-11 


3.92420373812 -1.0 x 10" 


-11 
















A^ 


/35 


-R35 




1000 


3.923939180351 


4.8 X 10-5 




2000 


3.924204540724 


8.0 X 10-^ 








4000 


3.924203740852 


-3.7 X 10-9 








6000 


3.924203738244 


7.9 X 10-11 








8000 


3.924203738121 


-4.4 X 10-12 








10000 


3.924203738135 


2.1 X 10-12 
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Table 2.2 



n 


tn 


n 


t-n 


n 


t-n 


1 





2 


4.1558656464 


3 


5.5491219562 


4 


8.1117756122 


5 


10.8952834492 


6 


12.0523651120 


7 


13.4542992617 


8 


14.9275108496 


9 


16.3036898019 


10 


17.7350625418 


11 


18.837088412 


12 


20.551890978 


13 


21.752187480 


14 


22.93715924 


15 


23.33859940 


16 


23.97767239 


17 


25.79365179 


18 


27.1212236 


19 


27.8904904 


20 


28.6462091 


21 


30.100668 


22 


30.884244 


23 


31.730116 


24 


32.248613 


25 


33.84677 


26 


34.08053 


27 


35.12990 


28 


36.04356 


29 


36.9637 


30 


38.2333 


31 


39.1512 


32 


39.7944 







Table 2.3 



A^ 


Rq 


/o R5 


h 


-Rio 


ho 


1000 


-3.04 -7.7 


X 10-2 -2.9558 1.4 x 10 


-2 -2.95033 


9.3 X 


10-4 


2000 


-2.85 3.1 


X 10-3 -2.9629 -3.5 x 10 


-3 -2.95651 


5.4 X 


10-4 


4000 


-2.95 -9.7 


X 10-3 -2.9588 7.3 x 10 


-4 -2.95642 


-8.7 X 


10-^ 


6000 


-2.95 6.3 


X 10-2 -2.9561 -1.4 X 10 


-3 -2.95659 


9.4 X 


10-^ 


8000 


-2.98 6.9 


X 10-2 -2.9558 -9.9 x 10 


-4 -2.95660 


4.8 X 


10-^ 


10000 


-2.97 1.8 


X 10-2 -2.9567 -2.3 x 10 


-4 -2.95659 


1.1 X 


10-^ 
















A^ 


-Rl5 


/l5 


-R20 


ho 




1000 


-2.956381 - 


-2.8 X 10-3 - 


-2.9591449 - 


4.1 X 10-4 




2000 


-2.956387 - 


-3.6 X 10-4 - 


-2.9568292 - 


9.6 X 10-^ 






4000 


-2.956613 - 


-2.3 X 10-^ - 


-2.9565854 


4.1 X 10-6 






6000 


-2.956592 - 


-3.8 X 10-*^ - 


-2.9565951 


2.5 X 10-6 






8000 


-2.956593 - 


-1.7x10-^ - 


-2.9565940 


6.7 X 10-^ 






10000 


-2.956593 - 


-2.3 X 10-^ - 


-2.9565932 


1.2 X 10-^ 


















A^ 


-R25 


/25 


-R30 


-^30 




1000 


-2.9584634 


3.5 X 10-3 


-2.95201141 


3.7 X 10-3 




2000 


-2.9566395 


1.9 X 10-4 


-2.95638520 


1.7 X 10-^ 






4000 


-2.9565937 


-7.5 X 10-6 


-2.95659865 


2.3 X 10-6 






6000 


-2.9565925 


1.3 X 10-^ 


-2.95659358 


-2.0 X 10-^ 






8000 


-2.9565932 


-1.8 X 10-s 


-2.95659347 


4.3 X 10-s 






10000 


-2.9565934 


-1.0 X 10-^ 


-2.95659342 


1.3 X 10-s 
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N 

1000 
2000 
4000 
6000 
8000 
10000 



-R35 

-2.951829903 
-2.956587888 
-2.956589949 
-2.956593420 
-2.956593344 
-2.956593405 



/35 

1.3 X 10-4 
-1.9 X lO"'* 

3.4 X 10"^ 
1.6 X 10"^ 
2.0 X 10"^ 

-5.0 X 10-9 



Table 2.4 



n 


Un 


n 


Un 


n 


Un 


n 


Un 


1 
5 
9 
13 


2.3864500 
8.841633 

13.76869 

18.0078 


2 
6 

10 
14 


4.3752457 
10.605890 
15.2075 
18.874 


3 

7 
11 


6.0435487 
11.437474 
15.6182 


4 
8 
12 


7.571907 
12.76622 
16.9663 









Table 3.1 








A^ 


Rq 


/o 


R5 


h 


Rio 


ho 




1000 


-1.686 


8.7 X 10- 


-2 -1.77582 


1.6 X 10-2 


-1.782116715 1.7x10" 


-3 


2000 


-1.835 


-2.9 X 10" 


-2 -1.78432 


-1.4 X 10-3 


-1.783663963 -1.9 x 10" 


-5 


4000 


-1.758 


-6.5 X 10" 


-3 -1.78364 


-1.6 X 10-4 


-1.783642384 -1.0 x 10" 


-6 


6000 


-1.820 


-5.2 X 10" 


-3 -1.78366 


2.0 X 10-^ 


-1.783642948 8.6 x 10" 


-8 


8000 


-1.760 


-1.1 X 10" 


"2 -1.78364 


-2.0 X 10-^ 


-1.783642816 -2.7 x 10" 


-8 


10000 


-1.786 


8.2 X 10" 


-3 -1.78363 


1.1 X 10-^ 


-1.783642826 6.5 x 10" 


-9 


















A^ 


R 


L5 


/l5 


-R20 




ho 




1000 


-1.783149679719 


9.4 X 10"^ 


-1.7835293644621857 


-2.4 X 10-^ 




2000 


-1.783644577360 


-3.2 X 10-8 


-1.7836428890128260 


8.1 X 10-9 




4000 


-1.783642829247 


-2.2 X 10-*^ 


-1.7836428272010046 


-1.3 X 10-11 




6000 


-1.783642827253 


8.0 X 10"^^ 


-1.7836428271539760 


4.5 X 10-13 




8000 


-1.783642827162 


-6.1 X 10-^2 


-1.7836428271544534 


-7.0 X 10-1^ 




10000 


-1.783642827150 


-1.4 X 10-^2 


-1.7836428271544125 


-8.3 X 10-1^ 
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N 


-R25 


/25 






1000 


-1.7836271941181258672 - 


1.3 X 10-^ 




2000 


-1.7836428277732280167 


1.1 X 10-9 






4000 


-1.7836428271546296417 


1.4 X 10-1^ 






6000 


- 1 . 7836428271544148349 


9.0 X 10-1*^ 






8000 


-1.7836428271544160390 


2.0 X 10-is 






10000 


-1.7836428271544160184 - 


■6.2 X 10-19 














A^ 


-R30 


-^30 




1000 


-1.7836430434865235424107 


-3.5 X 10-6 




2000 


-1.7836428270926205946250 


5.0 X 10-11 






4000 


-1.7836428271544153839298 


1.1 X 10-1^ 






6000 


-1.7836428271544160178686 


-1.6 X 10-is 






8000 


-1.7836428271544160181358 


2.7 X 10-20 






10000 


-1.7836428271544160181708 


-3.8 X 10-22 














A^ 


-R35 


/35 




1000 


-1.7836435960725853092305493 -2.3 x 10-^ 




2000 


-1.7836428271502571300252154 -3.7 x lO-^^ 






4000 


-1.7836428271544160011102849 4.2 x 10-^^ 






6000 


-1.7836428271544160181784305 -8.5 x lO-^i 






8000 


-1.7836428271544160181689502 1.0 x 10-^^ 






10000 


-1.7836428271544160181690199 2.9 x 10-^^ 








Table 3.2 




N 


Rq Iq -R5 


-^5 -Rio 


/i( 


) 


1000 


3.43 -3.9 X 10-1 3.3563 


-3.5 X 10-1 3.2408605 


-3.5 X 


10-1 


2000 


3.23 8.9 X 10-2 3.1032 


6.7 X 10-3 3.0945490 


2.3 X 


10-4 


4000 


3.10 2.6 X 10-2 3.0906 


-3.1 X 10-4 3.0913811 


2.3 X 


10-6 


6000 


3.12 4.1 X 10-2 3.0912 


3.7 X 10-4 3.0914333 


2.3 X 


10-6 


8000 


3.07 -4.2 X 10-3 3.0914 


-1.7 X 10-4 3.0914340 


-3.9 X 


10-^ 


10000 


3.12 -1.8 X 10-2 3.0914 


-1.3 X 10-4 3.0914342 


-1.6 X 


10-^ 












A^ 


-R15 A5 


-R20 ho 




1000 


3.0854245832 -3.1 x lO-i 


2.923318726655 -2.4 x 


10-1 




2000 


3.0919697213 -3.7 x 10-^ 


3.091486059342 -6.5 x 


10-5 




4000 


3.0914327625 -3.9 x 10-^ 


3.091434224098 4.0 x 


10-8 




6000 


3.0914342708 6.1 x 10-^ 


3.091434237097 5.7 x 


10-10 




8000 


3.0914342301 -3.7 x 10-^ 


3.091434236714 -1.0 x 


10-11 




10000 


3.0914342356 -3.6 x 10-i° 


3.091434236735 -1.3 x 


10-12 
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R 



N 

1000 
2000 
4000 
6000 
8000 
10000 3.09143423673865 



25 

2.80762049069122 
3.09143198592104 
3.09143423720919 
3.09143423674457 
3.09143423673861 



-^25 

-1.1 X 10- 

-2.2 X 10- 

2.4 X 10- 

-5.1 X 10 

2.3 X 10 

1.0 X 10 



-1 

-5 

-9 

-12 

-13 

-14 



^30 

2.78668118045678635 
3.09143105396801472 
3.09143423678556380 
3.09143423673870285 
3.09143423673865144 
3.09143423673865098 



ho 

7.4 X 10-2 
-3.5 X 10-6 

3.9 X 10-11 
-1.7 X 10-13 

1.9 X 10-1^ 

6.2 X 10-1^ 



A^ 



R 



35 



^35 



1000 2.91398635132232829470 

2000 3.09143341359018088974 

4000 3.09143423673979921383 

6000 3.09143423673865010310 

8000 3.09143423673865095777 

10000 3.09143423673865095240 



2.7 X 10-1 
1.4 X 10-^ 
-2.4 X 10-12 
-1.1 X 10-1^ 
-1.9 X 10-18 
-1.3 X 10-20 



Table 3.3 



n 


'■n 


n 


tn 


n 


'■n 


1 


8.2720409199 


2 


11.3959869930 


3 


14.8616932015 


4 


17.1783243050 


5 


19.2124566315 


6 


20.8274294554 


7 


23.4659374198 


8 


25.2726883522 


9 


27.0035774491 


10 


28.1569222690 


11 


30.2145623343 


12 


31.6193141164 


13 


33.7856279775 


14 


34.9435854723 


15 


36.5559515067 


16 


37.6356026748 


17 


39.1608229256 


18 


40.6589300308 


19 


42.8804581030 


20 


43.2736012304 


21 


44.9765395474 


22 


46.4176568046 


23 


47.2517710599 


24 


48.7821808287 


25 


50.3519022325 


26 


51.5688981695 


27 


53.1356287828 


28 


54.0717837181 


29 


55.0990003336 


30 


56.4089955139 


31 


57.5391214415 


32 


59.1986375433 


33 


60.1739007171 


34 


61.6441827270 


35 


62.8146545420 


36 


63.4247884022 


37 


65.1023702197 


38 


66.0180646898 


39 


66.8050237006 


40 


68.6802278238 


41 


69.8132058342 


42 


70.7502185552 


43 


71.9861530156 


44 


72.7927328082 


45 


74.1137296216 


46 


74.8761895173 


47 


76.2796967025 


48 


77.4608764665 


49 


78.7319975717 


50 


79.5372511477 


51 


80.8499015926 


52 


81.9286308045 


53 


82.6995529553 


54 


83.4681179192 


55 


85.2402769759 


56 


85.6802121224 


57 


87.2830188249 


58 


88.5094955323 


59 


89.2377130355 


60 


90.0534073382 


61 


91.4472572430 


62 


92.0496894589 


63 


93.3566370961 


64 


94.2221147468 


65 


95.3044565474 


66 


96.6527715250 


67 


97.7264314003 


68 


98.4244540180 


69 


99.4730638315 
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Table 3.4 



n 


Un 


n 


Un 


n 


Un 


1 


12.9399446108 


2 


15.1248640287 


3 


17.2775088490 


4 


21.9119654118 


5 


23.7124474310 


6 


27.6868648494 


7 


29.1470584255 


8 


31.1315265360 


9 


31.9862854000 


10 


33.6323734231 


11 


35.7361264638 


12 


38.1008875317 


13 


39.9548075690 


14 


41.3629251312 


15 


43.0030848131 


16 


43.7924301232 


17 


49.3874980802 


18 


50.3892911690 


19 


51.9883497256 


20 


53.3610715851 


21 


55.5058736308 


22 


57.1306190068 


23 


58.5145765119 


24 


59.2810504632 


25 


60.8114113807 


26 


61.7177742037 


27 


62.3299217969 


28 


65.1200148215 


29 


66.3871768599 


30 


67.7658152255 


31 


68.5636482576 


32 


70.0994795387 


33 


71.9139076205 


34 


73.4598285562 


35 


74.4219698604 


36 


75.9259426071 


37 


76.7219513797 


38 


80.2370604179 


39 


80.9795116625 


40 


82.2141987387 


41 


84.1686266809 


42 


85.4558525934 


43 


86.4596989555 


44 


87.5017176801 


45 


88.7989284271 


46 


90.8930274655 


47 


91.6094970880 


48 


93.0648834307 


49 


93.8455295242 


50 


94.4090447654 


51 


95.8637476000 


52 


96.5025102751 


53 


97.8219970593 


54 


98.9086789539 



Table 3.5 



n 


Pn 


n 


Pn 


1 


3.2308208282 + 8.9496290911 i 


2 


3.0144204971 + 19.1670355895 z 


3 


3.1880664988 + 26.3033849287 i 


4 


3.1549639910 + 36.6242398231 i 


5 


2.7150409653 + 45.1719799932 i 


6 


2.4938210677 + 47.5816502442 i 


7 


3.3624175212 + 54.4320525502 i 


8 


2.9077749773 + 64.2513434784 z 


9 


3.1556119321 + 71.9344926377 i 


10 


2.4066868777 + 78.1144688947 z 


11 


2.9102154501 + 82.3890698796 i 


12 


2.8870784016 + 89.7875787849 z 


13 


3.3672596002 + 99.9194124003 i 


14 


2.7073645379 + 107.1592978688 z 


15 


2.7721770492 + 110.2613188689 i 


16 


3.1645227049 + 117.3896482956 z 


17 


2.5547542302 + 126.5001914198 i 


18 


2.8669588498 + 128.2070045099 z 


19 


3.1444059195 + 135.3354942155 i 


20 


3.3550056102 + 145.3938901873 z 
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Table 3.6 



n 


«n 


n 


«n 


n 


«n 


1 

4 

7 

10 

13 

16 

19 

22 

25 

28 

31 

34 

37 

40 

43 

46 

49 

52 

55 

58 

61 

64 

67 

70 

73 

76 


5.6185671952 
13.7962357292 
19.4575970308 
24.1962387103 
29.1281915898 
33.0414393853 
36.6630191145 
40.6779144658 
44.5137568744 
47.8221666146 
51.5155718913 
54.4042145877 
59.7418610167 
63.2789136582 
66.3286609040 
71.4492522675 
74.5424808567 
77.0682809138 
80.0330330946 
83.0719604785 
85.7778751490 
88.7059787721 
91.4583733416 
94.5091774704 
96.9317027531 
99.8881597950 


2 

5 

8 

11 

14 

17 

20 

23 

26 

29 

32 

35 

38 

41 

44 

47 

50 

53 

56 

59 

62 

65 

68 

71 

74 


9.3587692608 
16.0894994874 
21.8594962784 
26.0225517432 
29.9334498107 
33.9453293541 
38.2010379041 
41.3507131813 
45.2712575072 
49.0105055856 
51.6631984921 
55.5750957955 
61.2006909669 
64.1010371646 
67.5430514279 
72.3027663172 
75.2004527728 
78.2931224271 
81.1315223512 
84.1964169499 
86.9407183317 
89.9368381341 
92.7393858316 
95.4769845208 
97.8372552255 


3 

6 

9 

12 

15 

18 

21 

24 

27 

30 

33 

36 

39 

42 

45 

48 

51 

54 

57 

60 

63 

66 

69 

72 

75 


12.0264936925 
17.7847280762 
22.5758373316 
27.3617234087 
31.7848563053 
35.7534640159 
39.7144586843 
43.0237415354 
46.4744214908 
49.9881789218 
53.5375246718 
58.8270009777 
61.7440856002 
64.8912913186 
70.1770012608 
73.1469069994 
76.3094730233 
79.2581319133 
81.9704731391 
85.2011952740 
87.4479355411 
90.6181312264 
93.0881396596 
96.1116129723 
99.1010053953 


n 


Un 


n 


Un 


n 


Un 


1 

4 

7 

10 

13 

16 

19 

22 

25 

28 

31 

34 



29.3670678246 
38.5418813017 
46.8484465576 
54.2746979473 
58.8640673277 
65.1978869599 
75.0463737748 
80.0462305996 
86.6795992346 
92.9081333908 
97.7100184923 


2 

5 

8 

11 

14 

17 

20 

23 

26 

29 

32 

35 


24.4022873037 
33.5747954436 
40.0447318001 
50.0699839799 
55.4033486176 
60.9985910184 
71.4691922483 
76.6579317410 
83.8815457436 
88.7015447955 
94.5315054431 
99.3454122963 


3 

6 

9 

12 

15 

18 

21 

24 

27 

30 

33 


26.6418851276 
35.2863556538 
44.8539372315 
52.1256566323 
56.5807632829 
63.8456398387 
72.8978373808 
78.6419254744 
85.2027989008 
90.5599843766 
96.4777165784 
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Table 3.7 



n 


Pt 


n 


Pt 


1 


3.3591319232 + 57.3250633340 z 


2 


3.5001503209 + 68.5679322965 i 


n 


Pn 


n 


Pn 


1 

3 

5 

7 

9 

11 

13 

15 

17 

19 

21 

23 


5.7692647648 + 8.9956889852 z 
5.3846794177 + 18.2757545274 i 
5.3861276711 + 26.6587658619 z 
5.7862146670 + 36.6524415925 z 
5.2171816339 + 44.9332286844 z 
5.6527005509 + 54.2712676115 z 
5.2007749392 + 63.8919812818 z 
3.5790185110 + 69.3185823254 z 
4.6599017691 + 76.8392451729 z 
3.3026675429 + 82.2734272886 z 
4.9626158230 + 90.6464501774 i 
5.8412001878 + 99.7795574645 z 


2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

22 


4.8476735625 + 14.0858508094 i 
4.2067408135 + 20.6821222248 z 
4.5402410961 + 31.5997171480 z 
4.3281324452 + 41.9810464020 z 
4.7664043254 + 48.8240853600 z 
4.6702375442 + 60.1227186672 z 
3.7792244502 + 66.9023711945 z 
5.7686816361 + 72.2279984677 z 
5.3421407237 + 81.9810448650 z 
3.9986981918 + 88.8368300771 i 
4.7401362618 + 94.8547536615 z 



Table 4.1a 



A^ 


-Ro -^0 


Ri h 


-R2 I2 




1000 


2.22 6.3 X 10 


-3 2.1777 -1.1 X 10-2 


2.173977 2.6 x 10-4 




5000 


2.14 -2.2 X 10 


-2 2.1700 3.7 X 10-3 


2.173544 -1.2 x 10-4 




10000 


2.18 -1.3x10 


-2 2.1753 1.1 X 10-3 


2.173833 -1.0 X 10-4 




30000 


2.19 -2.2 X 10 


-2 2.1729 -6.9 X 10-4 


2.173759 -1.0 X 10-^ 




100000 


2.17 1.2 X 10 


-2 2.1735 1.0 X 10-4 


2.173747 -4.7 x 10-^ 














A^ 


R3 


/3 Ra L 


I -R5 -^5 




1000 


2.167862 -6.2 


X 10-3 2.20484 -5.4 x 


10-2 2.49186 7.1 X 10" 


-2 


5000 


2.173880 -6.8 


X 10-4 2.17935 -1.6 x 


10-4 2.18302 2.7 x 10" 


-2 


10000 


2.173937 1.8 


X 10-4 2.17267 1.6 x 


10-3 2.16380 -2.8 x 10" 


-3 


30000 


2.173700 -2.3 


X 10-6 2.17367 -3.7 x 


10-4 2.17551 -9.3 X 10" 


-4 


100000 


2.173745 -5.9 


X 10-6 2.17378 -5.1 x 


10-^ 2.17406 9.7 x 10" 


-5 
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N 

1000 

5000 

10000 

30000 

100000 



Rq 

2.352 
2.082 
2.170 
2.179 
2.173 



1.2 

7.3 X 10-2 

-4.1 X 10-2 

5.3 X 10-3 

1.3 X 10-3 



R7 
-1.442 
1.835 
2.293 
2.165 
2.169 



2.2 

-1.7 X 10-1 

-6.2 X 10-2 

2.4 X 10-2 

1.8 X 10-3 



Rs 
-7.975 
2.195 
2.480 
2.101 
2.164 



-6.1 

-1.0 

2.2 X 10-1 

1.0 X 10-2 

-7.2 X 10-3 



Table 4.1b 



A^ 



Rq 



'0 



R 



h 



Ri 



h 

2.6 X 10-2 
-8.6 X 10-4 

6.8 X 10-4 

1.7 X 10-4 

4.9 X 10-5 



1000 

5000 

10000 

30000 

100000 



N 



3.263 
3.305 
3.245 
3.263 
3.276 



8.3 X 10-2 
-1.1 X 10-2 

6.4 X 10-3 
2.9 X 10-2 

-8.2 X 10-3 



3.3048 
3.2708 
3.2612 
3.2697 
3.2665 



6.1 X 10-2 

-1.8 X 10-2 

-4.8 X 10-3 

9.0 X 10-4 

5.8 X 10-4 



3.28176 
3.26879 
3.26487 
3.26616 
3.26610 



R:^ 



Ra 



R, 



1000 

5000 

10000 

30000 

100000 



3.38022 
3.26878 
3.26486 
3.26672 
3.26609 



2.5 X 10-2 

6.2 X 10-3 

-2.7 X 10-3 

-1.6 X 10-4 

6.8 X 10-5 



3.524 
3.080 
3.329 
3.271 
3.264 



2.2 

7.2 X 10-2 
-3.1 X 10-2 

1.2 X 10-2 
-1.3 X 10-3 



-33.06 
1.75 
3.88 
3.06 
3.28 



8.3 
-2.9 
1.0 
1.0 X 10-1 
-3.1 X 10-2 



Table 4.2 



n 


^n 


n 


tn 


n 


tn 


n 


^n 


1 


10.2562 


2 


13.6866 


3 


15.9599 


4 


17.038 


5 


19.026 


6 


20.017 


7 


22.472 


8 


23.745 


9 


25.351 


10 


26.229 


11 


27.561 


12 


28.847 


13 


29.986 














-1 


-3.12740 


-2 


-6.5577 


-3 


-8.8310 


-4 


-9.9092 


-5 


-11.8976 


-6 


-12.888 


-7 


-15.343 


-8 


-16.616 


-9 


-18.222 


-10 


-19.101 


-11 


-20.433 


-12 


-21.718 


-13 


-22.857 


-14 


-24.859 


-15 


-25.892 


-16 


-26.850 


-17 


-28.418 


-18 


-28.927 
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Table 4.3 



k 


n 


T 


N{T) 


k 


n 


T 


N{T) 


Q(v^) 


1 


100 


84 


Q(v^) 


-1 


100 


93 


Q(v^) 


2 


100 


82 


Q(v^) 


-2 


100 


96 


Q(v^) 


3 


100 


81 


Q(v^) 


-3 


100 


98 


Q(v^) 


4 


100 


80 


Q(v^) 


-4 


100 


100 


Q(v^) 


5 


100 


78 


Q(v^) 


-5 


100 


101 


Q(v^) 


10 


100 


79 


Q(v^) 


-10 


100 


108 


Q(v^) 


1 


100 


75 


Q(^) 


-1 


100 


88 


Q(v^) 


2 


100 


72 


Q(v^) 


-2 


100 


92 


Q(v^) 


3 


100 


71 


Q(v^) 


-3 


100 


95 


Q(v^) 


4 


100 


71 


Q(v^) 


-4 


100 


97 


Q(v^) 


5 


100 


71 


Q(v^) 


-5 


100 


99 


Q(v^) 


10 


100 


93 


Q(v^) 


-10 


100 


106 


Q(V19) 


1 


50 


51 


Q(V19) 


-1 


50 


52 


Q(Vl9) 


2 


50 


50 


Q(v/19) 


-2 


50 


53 


Q(V29) 


1 


100 


107 


Q(v/29) 


-1 


100 


112 


Q(V29) 


2 


100 


105 


Q(v/29) 


-2 


100 


114 


Q(v/29) 


3 


100 


104 


Q(v/29) 


-3 


100 


116 


Q(V29) 


4 


100 


102 


Q(v/29) 


-4 


100 


117 


Q(v/29) 


5 


100 


102 


Q(v/29) 


-5 


100 


118 


Q(V29) 


10 


100 


99 


Q(V29) 


-10 


100 


123 


Q(V3i) 


1 


40 


41 


Q(v/31) 


-1 


40 


42 


Q(V67) 


1 


30 


32 


Q(V67) 


-1 


30 


33 



34 









Table 5.1 






A^ 


Rq 


/o 


-R2 h 


R4 


h 


1000 


3.09 


5.3 X 10-1 


3.334 4.5 X 10 


-1 3.483 


4.9 X 10-1 


5000 


3.31 - 


-5.5 X 10-2 


3.317 -2.4 X 10 


-2 3.335 


-3.2 X 10-3 


10000 


3.42 


1.5 X 10-1 


3.417 -7.8 X 10 


-3 3.389 


-1.3 X 10-2 


30000 


3.45 - 


-1.4 X 10-1 


3.385 -2.3 X 10 


-2 3.380 


-6.2 X 10-3 


100000 


3.25 - 


-4.8 X 10-2 


3.375 1.6 X 10 


-3 3.383 


1.2 X 10-3 




N 


Rq 


h 


Rs 


/8 


-Rio -^10 


1000 


3.9634 


6.4 X 10" 


■1 4.03999 6.2 


X 10-1 4.792707 5.2 x lO-i 


5000 


3.3419 


3.6 X 10" 


"2 3.38817 4.1 


X 10-2 3.400866 6.3 x IO-2 


10000 


3.3825 


-1.0 X 10" 


"2 3.37107 -4.5 


X 10-3 3.377550 4.8 x 10-^ 


30000 


3.3803 


1.9 X 10" 


"4 3.38367 1.4 


X 10-3 3.383935 3.6 x 10-^ 


100000 


3.3839 


1.2 X 10" 


"4 3.38369 -8.5 


X 10-^ 3.383657 -8.2 x 10-^ 








Table 5.2 







n 


Un 


n 


Un 


n 


Un 


n 


Un 


n 


Un 


1 
6 


2.79373 
7.46 


2 

7 


4.0887 
7.90 


3 

8 


5.362 

8.8 


4 
9 


5.887 
9.6 


5 


7.03 



Table 6.1 



N 










N 


(5) 


N 




5 


0.96 


0.98 


0.93 


0.92 





1.03 





0.98 


10 


0.95 


1.05 


1.12 


0.84 


2 


0.82 


1 


0.99 


15 


1.07 


1.06 


0.85 


1.03 


4 


0.81 


2 


1.01 


20 


1.07 


1.09 


1.23 


1.07 


6 


1.19 


3 


0.99 


25 


1.06 


0.90 


1.25 


0.87 


8 


1.39 


4 


1.01 


30 


0.88 


0.88 


0.80 


1.22 


10 


0.67 


5 


0.94 
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Figure 7.1 
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